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. We compute the next-to-leading order (NLO) QCD corrections to the spin-independent and 

O ' spin-dependent cross sections for the production of a single-hadron accompanied by an opposite 

O ■ jet in hadronic collisions. This process is being studied experimentally at RHIC, providing a 

, new tool to unveil the polarized gluon distribution Ag. We perform a detailed analysis of the 

^ ■ phenomenological impact of the observable at NLO accuracy and show that the preliminary data 



by the STAR collaboration confirms the idea of a small gluon polarization in the 0.05 ^ x < 0.3 
range. 



PACS numbers: 13.88. +e, 12.38.Bx, 13.87.Fh 



April 2009 



2 



I. INTRODUCTION 

In the past two decades, measurements of the spin asymmetries {N = p, n, d) in longitu- 
dinally polarized deep-inelastic scattering (DIS) [l| have provided new information on the spin 
structure of the nucleon. One of the most surprising results is that only a small fraction of the 
spin of the proton can be attributed to the spin of the quarks. The main goal of the spin program, 
besides obtaining the partonic share to the total spin of the nucleon, is the extraction of the full 
set of the x— dependent polarized quark (Ag = — q^) and gluon (A(7 = — g^) densities of 
the nucleon. Many phenomenological analyses [2! demonstrate, however, that available DIS data 
alone is not sufficient for this purpose. This is true in particular for ^g{x^ Q"^) since it contributes 
to DIS in leading-order (LO) only via the dependence of gi (or Ai) which could not yet be 
accurately studied experimentally. As a result of this, it turns out that the x shape of A^f seems 
to be hardly constrained by the DIS data. 

The precise extraction of Ag thus remains one of the most interesting challenges for spin physics 
experiments. The RHIC collider at BNL, running in a proton-proton mode with longitudinally 
polarized beams provides the ideal tool for that purpose. The observables measured so far include 



single-pion production at center-of-mass energy = 62 [31] and 200 GeV and jet production 
at ^/s = 200 GeV j^. Unlike DIS, those processes have a direct gluonic contribution already at 
the lowest order. 

A next-to-leading order (NLO) global analysis that includes all available data from inclusive 
and semi-inclusive polarized deep-inelastic scattering, as well as from polarized proton-proton 
scattering at RHIC has been recently performed j^. The main outcome of the analysis is the 
indication of a rather small gluon polarization in the nucleon over the limited region of momentum 
fraction 0.05 < x < 0.2. 

Q 

Recently the STAR collaboration at RHIC has presented preliminary data [7| corresponding 
to the 2006 run on a less inclusive observable, involving the production of a charged hadron 
accompanied by a back-to-back jet. From the pure experimental point of view, counting with an 
opposite jet allows one to use it as a trigger for the hadron, reducing the bias in the selection. 
Furthermore, having a more exclusive observable, and particularly counting with the transverse 
momentum of both the hadron and the jet, permits one to perform a more detailed study to 
extract Ag. 

In order to make reliable quantitative predictions for a high-energy process, it is crucial to 
determine the NLO QCD corrections to the Born approximation. In general, in hadronic collisions. 
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cross sections computed at the lowest order in perturbation theory are severally affected by the 
dependence on the 'unphysical' factorization and renormalization scales, dependence that can be 
partially cured only by including the NLO corrections. Furthermore, the appearance of one extra 
final-state parton in the NLO from the 2 —> 3 real corrections allows one to improve the matching 
between the theoretical calculation and the realistic experimental conditions, particularly when 
jets are present. 

The calculation of the NLO QCD corrections to hadron+jet production by unpolarized and 
polarized hadrons is the purpose of this paper. Several modern versions of the subtraction method 
to calculate any infrared-safe quantity in unpolarized collisions are at present available in the 



literature j^, . The formalism of Ref . [9|] has been used in Ref . [10|] to construct a Monte Carlo 
code that can calculate any jet infrared-safe observable in hadron-hadron unpolarized collisions 
and generalized, in Refs. H, 12], to the polarized case. In isl, the method was extended to the 
case of singe-hadron inclusive observables. 

In the present paper, we apply the method of Refs. Il3l ] to the case of hadron-|-jet 



observables. As a result, we will present a customized code, with which it will be possible to 
calculate any infrared-safe quantity corresponding to one-hadron+jet production to NLO accuracy, 
for both polarized and unpolarized collisions. With the technique introduced in the code is 
suited to allow the inclusion of the observable in a global fit analysis. 

With the code at hand, we analyze in detail the phenomenological implications of the observ- 
able. The key point here is that, being a more exclusive measurement, it is possible to analyze the 
data in terms of a new set of variables and, also, to impose experimental cuts that can enhance the 
contribution of some partonic subprocesses over others. That is fundamental in order to increase 
the sensitivity on the spin-dependent gluon distribution in polarized collisions. 

This paper is organized as follows: in section 2 we describe the main ingredients of the calcula- 
tion and study the perturbative stability of the NLO results, by looking at the scale dependence 
and 'i^'-factors'. In section 3 we discuss some phenomenological aspects of hadron+jet production 
in hadronic collisions and in section 4 we present the results for the asymmetries at RHIC . Finally, 
section 5 contains the conclusions. 
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II. NLO CORRECTIONS AND PERTURBATIVE STABILITY 



The factorization theorem [1^ allows one to write the cross section for one-hadron production 
in hadronic collisions as 



daPP-h+jetx ^ jdxidx^dz f^'{x,,filj) ft {^2, fill) 



fuf2j 

X da-'^'-'^^^-'^^' {XiPi,X2P2,Ph/z, HfI, I^FF, fJ'R) D'}{z,fxlp) X S{pi) (1) 

where, Hi and H2 are the colliding particles with momentum pi and p2-, respectively, h is the 
outgoing hadron with momentum ph and the sum in Eq.([l]) runs over all possible initial and final 
partonic states. The parton distributions ff" are evaluated at the factorization scale ^pii the 



fragmentation functions at the scale ^ff [23| and the coupling constant, appearing in the pertur- 
bative expansion of the partonic cross section, at the renormalization scale yU^. The measurement 
function S{pi) accounts for possible experimental cuts applied to the cross section, and in this par- 
ticular case, for the definition of the jet in terms of the kinematics of the final-state partons. The 
analogous of Eq.([l]) for polarized cross section is obtained by replacing the parton distributions 
and the partonic cross section by its polarized expressions, A/^ ' and rfAo"-'^'^^^" , respectively. 
As usual, the (longitudinally polarized) asymmetry is defined by the ratio between the polarized 
and unpolarized cross sections 

In order to evaluate the QCD corrections to the process we rely on the version of the subtraction 



method introduced and extensively discussed in Refs. [9I, |10| , and in Ref . [Uj . We refer the reader 



n 



to those references for the details. The implementation is performed in a MonteCarlo like code 



profiting from the one available for the computation of single-hadron production in Ref. 13|. Since 



the calculation in [l3j provides the full kinematics for the final-state partons it is possible, after 
some modifications, to build the jet kinematics in terms of them and obtain the required cross 
section. It is worth noticing that the same code computes both unpolarized and polarized cross 
sections, since the formal structure of the corrections is exactly the same. 

In this work, we concentrate on the phenomenology of pion production accompanied by a back- 
to-back jet for the kinematics of the STAR experiment at RHIC with a center-of-mass energy of 
= 200 GeV. Unless otherwise stated, we require the pion transverse momenta to be larger 
than 2 GeV and the one for the jet to 10 GeV < p^^* < 25 GeV. The rapidities of the pion and 
the jets are limited to the range |?7| < 1. They have to be separated in the azimuthal angle 
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by A0 = — > 2 to ensure the pion and jets are produced from 'opposite-side' partons. 
Finally, the jets are defined according to the cone algorithm with R=0.7. 
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FIG. 1: Unpolarized (solid) and polarized (dashes) NLO K factors. The choice of the factorization and 
renormalization scales corresponds to fip = fJ-R = (^Pt Pt^^ Z^- 

The size of radiative QCD corrections to a given hadronic process is often displayed in terms 
of a i^-factor which represents the ratio of the NLO over LO results. In the calculation of the 
numerator of K, one obviously has to use NLO-evolved parton densities. As far as the denominator 
is concerned, a natural definition requires the use of LO-evolved parton densities. In the polarized 
case, a problem arises for such a definition: since the polarized pdfs are not as well constrained 
as the unpolarized ones, it might happen that quite different results for the A/'s can emerge 
when the fit is performed at LO or at NLO. This is particularly enhanced by the fact that the 
polarized pdfs have nodes. Therefore, the /^-factor for a given process, defined using LO parton 
densities in the denominator could be largely affected by the fact that some polarized densities 
are at present not well constrained. In order to avoid this problem, we define the /T-factor as the 
ratio between the NLO and the 'Born' cross section, where the latest corresponds to the use of 
NLO-evolved parton densities (and two-loop expression for as) when evaluating the lowest-order 
partonic cross sections in the denominator. Nevertheless, it is important to remember that the 
i^-factor is not a physical quantity and just provides a number to 'quantify' the effect of the 
higher-order corrections. 

In the unpolarized case, we use the MRST2002 parton distributions Differences of the 
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order of percent in the cross section are observed when more recent distributions are considered. 
Nevertheless, since the DSSV {g] distribution we use in the polarized case set is obtained from a 
global analysis that relies on the unpolarized MRST2002 set as a reference, we will restrict the 
analysis to the MRST2002 densities. For the fragmentation functions we rely on the DSS [l^ 
set that provides full flavor and charge separation at NLO. Similar results for pions are obtained 
with the latest AKK set [l^]. Considering that the process depends on two different hard scales, 
the transverse momentum of the hadron and the one of the jet, we define the 'default' scale to 
be the average of both of them, and, unless otherwise stated we set all the factorization and 
renormahzation scales to fip = fJ'R = {pt ~^ Pt^) /2- 

Fig. [1] shows the result for the unpolarized (solid) and polarized (dashes) 'i^— factors' computed 
at the default scales in terms of Px/Pt^ i ratio between the transverse momentum of the pion 
and the jet 3]- The relevance of this dimensionless ratio will become clear in the next section. 
As can be observed, the NLO corrections are sizable in the unpolarized case, ranging from 50% 
to more than 100% of the Born result. The corrections are smaller in the polarized case, as 
usual, resulting in an important decrease in the corresponding spin asymmetry. In that sense 



the situation is quite similar to the one found for single-hadron production 
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18| . In some 
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FIG. 2: LO and NLO unpolarized (left) and polarized (right) cross sections at different factorization and 
renormahzation scales = = ^pj, + I^x^^ results were scaled by a factor of 10""'^ /( 

better presentation. 



or 



extreme kinematical regimes, the QCD corrections tend to be dominated by the contributions 
arising from soft-gluon emission that need to be resummed to all orders in the coupling constant 
as to allow for a quantitative study. It has been shown [l9] that for inclusive single-hadron 
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production the resummation of the dominant terms is required in the case of fixed-target energy 
experiments, where the 'i^-factors' largely exceed those found here, but not for a collider running at 
y/s = 200 GeV. Particularly, the effect of the resummation over the corresponding asymmetries is 



rather small 20|. Considering that the dominant soft contributions for the hadron+jet observable 
originate from the same Sudakov form factors as in the inclusive case, we believe those effects can 
also be neglected in a first approach here. 

A reliable error estimate on our NLO results requires some knowledge on the size of the un- 
calculated higher-order terms. The best we can do, before higher-order terms are computed, is 
to study the dependence of the full NLO results on the renormalization and factorization scales. 
Although physical observables are obviously independent of the scales, theoretical predictions do 
have such a dependence, arising from the truncation of the perturbative expansion at a fixed 
order in the coupling constant a^. A large dependence on the scales, therefore, implies a large 
theoretical uncertainty. In order to show how the scale dependence is substantially reduced once 
the next-to-leading order corrections are included we will compare to the Born result. For the 
sake of presentation we set all the scales to be equal, and vary them by a factor of 2 up and down 
with respect to the default choice, i.e, fip = fiR = fi (p^ + p^*^*) /2 with /i = 1/2, 1, 2. Fig. [2]plots 
the corresponding scale dependence of the unpolarized (left) and polarized (right) cross sections, 
were we observe a considerable reduction when the NLO corrections are included. Nevertheless, 
it is worth noticing that the scale dependence is still rather large at NLO for the unpolarized 
cross section, of the order of ±20% or more (compared to about ±80% at the Born level). The 
scale dependence is much smaller in the polarized case, even reaching the stage in which at some 
kinematics the NLO cross section evaluated at /i = 2 is larger than the one at /i = 1/2, opposite to 
the LO expectation. Since the uncertainty in the unpolarized cross section directly contributes to 
the one for the asymmetry, one might consider the convenience of using directly Aa, instead of the 
asymmetry, to extract the polarized parton distributions with a considerably better theoretical 
accuracy. 

III. PHENOMENOLOGY 

As discussed in the Introduction, counting with the jet kinematics allows one to impose cuts 
that enhance the relevance of some kinematical region in the momentum fractions Xi y an d z. 



That can be observed in Fig. [3], where we plot the average value of the momentum (x) 2J] and 



hadronization (z) fractions for unpolarized single-vr"'' production (solid) and the corresponding 
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one for 7r"'"+jet (dashes), in both cases in terms of the transverse momentum of the pion. While 
the averages are relatively similar for x, the change is quite clear for the fragmentation fraction 
z. Whereas for single-hadron production (z) ~ 0.6 — 0.7 remains almost constant over the full 
kinematical range, it shows a large variation, from 0.2 to 0.8, when the opposite jet is required. 
The rise of (z) with the transverse momentum of the pion can be easily understood from simple 



jet 



0.0 L 



5 10 15 

7T 

Pt 



FIG. 3: Average of the partonic momentum fractions x and z for single pion production (solid) and pion 
accompanied by a jet with pj:^ > 10 GeV (dashes). 

physical considerations. At the Born level, only two final-state partons are produced, with opposite 
transverse momentum. For that kinematics, the pion is produced by the fragmentation of one of 
the partons, while the jet is just formed by the other one. Therefore, the ratio between the 
transverse momentum of the pion and the one of the jet is exactly the hadronization fraction z. 
Once a jet cut is applied, selecting the transverse momentum of the pion is equivalent to selecting 
the fraction of momentum that is transferred from the parton in the hadronization process. It is 
worth noticing that at the Born level, counting with the jet and hadron kinematics allows one to 
fully reconstruct all the momentum fractions as 

z = ^ 

jet 

Ft 

xi = [j^T^ e^v>{'njet) + p't^ e^vi'nh)) / ^fs (3) 
X2 = [p't exp {-rijet)+ Pt^ exp ( -77^ ) ) / . 

While those relations are not valid at NLO, since one more parton can be radiated, still one can 
observe that there is a strong correlation between the 'real' momentum fractions (the arguments 
of the parton distributions and fragmentation functions in Eq.([l])) and those obtained from the 
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FIG. 4: Correlations between the 'measured' and 'real' x (left) and z (right). 



measured observables in Eq.Q. The correlations found for 7r+ production in unpolarized collisions 
are plotted in Fig. HI Considering a bin of size 0.05 (0.1) for x (z), we find that about 90% (60%) 
of the generated (weighted) events in the MonteCarlo implementation of the NLO corrections give 
the same value for the 'real' and 'measured' momentum fractions, at least in the kinematical range 
where their contribution to the cross section is dominant. 
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FIG. 5: NLO unpolarized cross section in terms of the 'measured' (solid) and 'real' (dashes) partonic 
momentum fractions X (left) and z (right) 



The situation is also visible when the cross section is plotted in terms of the same variables, as 
shown in Fig. [5l In the dominant range of 0.05 ^ x < 0.3, the agreement between the 'measured' 
cross section and the one obtained in terms of the 'real' momentum fraction is at the percent level. 
The use of the variables in Eq.Q can therefore allow for an accurate reconstruction of the initial 
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state momentum fractions Xi2- In the case of the z distributions, the differences between the 'real' 
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FIG. 6: Contribution to the cross section from the gg (solid), qg (dashes) and qq (dots) channels in terms 
of the transverse momentum of the pion (left) and the variable z = (right). 



and the 'measured' quantities can reach up to 10% — 15%. However, it still becomes quite useful 
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FIG. 7: Contribution to the cross section from the gg (solid), qg (dashes) and qq (dots) channels in terms 
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to plot the cross section in terms of it. This is mainly because, by selecting a range in z, one can 
enhance or decrease the contribution from some partonic channel due to the particular behavior of 
the fragmentation functions. This feature can be observed in Fig. [6l where we show the fractional 
contribution to the NLO unpolarized cross section from the gg, qg and qq initial state partonic 
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channels 25|]. If the cross section is analyzed in terms of the transverse momentum of the hadron, 
as it happens for the single-inclusive case, the pure gluonic channel contribution gg becomes only 
sizable at small values of pl^ and then decreases rapidly, making the cross section less sensitive 
on the gluonic content of the proton. The situation changes when the same results are studied in 
terms of the variable z. Here the gg channel shows a larger fractional contribution over the entire 
kinematical range at expenses of a suppression of the pure quark channels, that at most account 
for only 20% of the cross section, providing an ideal scenario to extract ^g in polarized collisions. 
Something similar occurs for the same observable plotted in terms of the variable x, as shown in 
Fig. [71 The analyses confirm that hadron+jet production in hadronic collisions, in terms of both 
dimensionless variables x and provides a clear source of information on the gluon distribution. 
In the next section we will look directly at the correspondent sensitivity on Af? in polarized pp 
collisions. 



IV. ASYMMETRIES AT RHIC AND SENSITIVITY ON A5 



In order to analyze the sensitivity of the process on the polarized gluon distribution, we will 
compute the NLO asymmetries with three different sets of spin- dependent densities: DSSV 



GRSV (standard) [2l|, and GS-C j22|. The corresponding NLO distributions at = 50 GeV^, 
a typical scale for this process, are shown in the left-hand side of Fig. [HI As can be observed, the 
expectations from the three sets are quite different. While the DSSV distribution corresponds to 
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FIG. 8: Polarized gluon density at = 50 GeV^ from different sets of polarized pdfs (left) and their 
ratios to the unpolarized distribution (right). 
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the best fit from tlie latest global analysis of all polarized data p], the GRSV set can be considered 
as an 'upper bound' for the allowed range of gluon densities. The GS-C set provides a distribution 
compatible with the requirement of a small gluon polarization in the range 0.05 ^ x < 0.3 but 
with a node in that region and a very different behavior at smaller x compared to the DSSV set. 
The right-hand side in Fig. [8] shows the ratio between the corresponding polarized distribution 
and the unpolarized MRST2002 set. In the quark sector, the dominant distributions are very 
similar among the different sets. Therefore, since the variations in the quark sector are much 
smaller than the ones for gluons, we can expect that any differences between predictions for the 
asymmetries that are found when using different polarized parton density sets are to be attributed 
to the sensitivity of the observable to A^f. 
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FIG. 9: Expected asymmetries for vr (left) and vr"^ (right) production at RHIC from different sets of 
polarized pdfs in terms of z. 



We start the presentation of the expected asymmetries by looking first at tt^ and ir^ production 
in terms of the variable z, as shown in Fig. [9l As expected, the asymmetries for the DSSV and 
GS-C distributions turn out to be small. When a set with a larger gluon distribution, like GRSV, is 
considered, the asymmetries increase to the l%-2% level. The asymmetries for positive pions show 
a stronger sensitivity on the polarized gluon distribution at large z, in line with similar findings for 
single pion production at large transverse momentum 13| . The differences and similitudes between 
negative and positive pion asymmetries can be easily understood: at small z, the 'favored' (like 
u —>■ 7r~^ and d —>■ 7i~) and 'unfavored' (as u ^ 7t~^ and d —>■ 7r~) fragmentation functions are 
rather similar and therefore the cross section is almost charge invariant in that range. On the 
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contrary, at larger ^, favored distributions overcome the unfavored ones. The asymmetries reflect 
the differences between the positive Am and negative Ad parton distributions in the polarized 
proton that contribute with a different weight to the cross section. In this regime, the AuAg 
channel becomes dominant for tt"*" production resulting in a larger sensitivity on the polarized 
gluon density. A similar analysis can be performed in terms of the variable x, as defined in 
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FIG. 10: Expected asymmetries for vr (left) and ir^ (right) production at RHIC from different sets of 
polarized pdfs in terms of x. 



Eq.([3]). Here, since z is integrated out, one can expect closer results for ir^ and tt^ production. 
This is observed in Fig. [101 where both asymmetries reflect the shape and order of the curves 
for the AG/G ratio plotted in the right-hand side of Fig. [HI Notice that the sensitivity on the 
gluon distribution for ir^ production is increased when the asymmetry is analyzed in terms of the 
variable x instead of z. 

The STAR collaboration at RHIC has recently presented preliminary data on vr+jet produc- 
tion. Even though the data has been obtained with very similar cuts to those used along this 
work, the experimental results can not be directly compared to the expected NLO asymmetries 
shown above because the data has not been corrected by the jet trigger efficiency. In order to 
allow for a comparison, we have recomputed the corresponding asymmetries by incorporating in 
the theoretical calculation the trigger efficiency parametrized as in Ref. [3] and modifying the 
experimental cuts accordingly. The only modification with respect to the previous cuts are those 
applying to the jet, which is required to have a transverse momentum in 9.5 GeV < p^"* < 25 GeV 
and rapidity in the range —0.7 < rj < 0.9 [2^. The result is presented in Fig. [TTJ The main effect 
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FIG. 11: Asymmetries measured by STAR at RHIC for 7r~ (left) and n'^ ( right) compared to the prediction 
from different sets of polarized pdfs. The theoretical predictions were corrected to account for the jet trigger 
efficiency. 

of the trigger is to enhance the contribution from large jp^^ with respect to the small p!^^ events 
and, therefore, increase the average (x) resulting in larger asymmetries. 

With the present experimental accuracy it is not yet possible to perform a precise extraction 
of the polarized gluon density from this observable. Nevertheless, the data can already rule out 
any possible scenario with a large gluon polarization in the range 0.05 ^ a; < 0.3. For that 
purpose we include in Fig. [11] the prediction from the set GRSV-max set, where the polarized 
gluon distribution is assumed to be equal to the unpolarized density at the very low intial scale of 
/i^ = 0.4 GeV^. That set completely overestimates the experimental data at small z. Therefore, 
in line with other measurements performed at RHIC, the preliminary data confirms the results 



from the global analysis in 



6| and points out to a small gluon polarization in the proton. 



V. CONCLUSIONS 



It is shown that the perturbative stability of the hadron+jet cross section improves considerably 
after including the NLO contributions. The corrections are found to be nontrivial: /T— factors are 
larger for the unpolarized cross section than for the polarized one, resulting in a reduction of the 
asymmetry at NLO. The possibility of looking at charged pions accompanied by a back-to-back 
jet is studied phenomenologically in detail, finding that the asymmetries for tt prodution, in terms 
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of both dimensionless variables x and are sensitive to the polarized gluon density in the range 
0.05 ^ x < 0.3. Furthermore, we show to NLO accuracy that the available data collected by 
the STAR collaboration at RHIC for this observable confirms the idea of a rather small gluon 
polarization. 
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